Quantum coherence engineering of mechanically motional states in optomechanical systems 
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We propose to synthesize arbitrary nonclassical motional states in optomechanical systems by using sideband 
excitations, which is an extensively used method in experiments, e.g., ground state cooling of the mechanical 
resonator in optomechanical systems. We first demonstrate how the Hamiltonian of the optomechanical systems 
can be reduced to that of trapped ions in the strong single-photon coupling regime when the photon blockade 
occurs, and then show a method to achieve target states. We further analyze the effect of the photon leakage 
on the fidelity of the expected state due to small anharmonicity. Moreover, we also discuss the experimental 
feasibility and provide operational parameters using the current technology. 
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Introduction. — Whether macroscopic mechanical res- 
onators behave quantum mechanics is a long-outstanding de- 
bate of the fundamental physics (see, e.g., reviews in Refs. jTl- 
Recent progresses of experimental techniques on, e.g., 
ground-state cooling and the fabrication of high-frequency 
mechanical resonators push forward the process to end this 
debate. In the existing literatures, several methods have been 
used to cool the mechanical resonators to their ground state 
in various types of the nano-structures, e.g., doubly clamped 
beams, singly clamped cantilevers, radial breathing mode of 
micro-toroids, and membranes. Moreover, the microwave- 
frequency mechanical resonator made by a thin film of alu- 
minium nitride has been experimentally demonstrated 
The potential applications of mechanical resonators in the 
quantum regime can be referred to, e.g., quantum informa- 
tion processing and quantum measurements such as sensitive 
detection of very weak forces. 

In quantum information processing, superpositions of 
quantum states are main resources. Many theoretical propos- 
als and experimental demonstrations have been presented to 
generate and manipulate quantum superposed states. For ex- 
ample, we have theoretically studied how to generate superpo- 
sitions of different Fock states for microwave photons jSfl, and 
later on experimentalists produced Fock states and also 
arbitrary superpositions of different Fock states by cou- 
pling a single-mode microwave cavity field to a superconduct- 
ing phase qubit. Similarly, particular non-classical phonon 
states of the vibration mode of trapped ions have been theo- 
retically studied (see, e.g., RefsJsUldl) and experimentally 
demonstrated (see, e.g., Refs. Ill ll 11211 '). The superpositions 
of macroscopically distinct quantum states have been theoret- 
ically studied in a mechanical resonator by coupling it to a 
charge qubit |, 1 3i1 . However, the generation of arbitrary non- 
classical motional states (hereafter, we call them as phonon 
states) in mechanical resonators is still an open question. 

Many approaches have been proposed to measure or con- 
trol mechanical resonators to the quantum ground-state by in- 
tegrating them with other quantum components, e.g, a single- 
mode cavity field or a superconducting qubit. Among these 
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FIG. 1: (Color online) (a) Schematic diagram for optomechanical 
systems with the radiation-pressure type interaction: the cavity can 
be in either the optical, or the microwave, or the radio-wave regime; 
and the mechanical resonator can be doubly clamped beams, singly 
clamped cantilevers, radial breathing mode of micro-toroids, and 
membranes, (b) Three different transition processes are presented 
by the black (carrier), red (fc-phonon red sideband excitation) and 
blue (fc-phonon blue sideband excitation) arrow lines. 



approaches, optomechanical systems, in which a macroscopic 
mechanical resonator is coupled to a single-mode cavity field 
through the radiation pressure, are extensively studied (see, 
e.g., review in Refs. H S). Recent experiments showed 
that the optomechanical systems are approaching the single- 
photon strong coupling regime lll4l4l7ll . In this regime, the 
photons transmitted through the cavity of the optomechanical 
system can be blockaded ITsl [l9ll and the single-photon op- 
tomechanical dynamics can be solved by using the stochastic 
Hamiltonian method |20]. 

In contrast to the proposals for generating the non-Gaussian 
quantum states (see, e.g., Refs. I'^l'-E?]), we here study a 
method to synthesize arbitrary nonclassical phonon states in 
strongly coupled optomechanical systems. Different from the 
recent proposal of measurement-based phonon state genera- 
tion 1 31], our method is deterministic and based on the side- 
band excitations, which have been theoretically fesi l26ll and 
experimentally studied to cool the mechanical resonator to its 
ground state [27-30]. We first demonstrate that the Hamil- 
tonian of the optomechanical system can be reduced to that 
of trapped ions ifs Uloll in the strong single-photon coupling 
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regime when the photon blockade occurs, and then show that 
the superpositions of phonon states can be generated using 
the similar approach as in trapped-ion systems. We further 
analyze the effect of the photon leakage on the fidelity of the 
target state. Moreover, we discuss the experimental possibil- 
ity and provide the operational parameters using the current 
technology. 

Theoretical model. — We study an optomechanical system, 
which can be either the membranes coupled to an optical cav- 
ity, or the optical cavity with one-end movable mirror, or the 
superconducting transmission line resonator coupled to a me- 
chanical beam. As schematically shown in Fig. \Tla), those 
systems have the radiation-pressure-type interaction, and the 
Hamiltonian of the optomechanical system driven by a classi- 
cal field can be written as 



H — hiOcd^a + hwmb^b + hga^ a{h'^ + b) 

+hn Lte-'('^«*+'^'') + h.c. 



(1) 



Here, (a) is the creation (annihilation) operator of the cavity 
field with the frequency lJc, and b^ (6) is the creation (annihila- 
tion) operator of the mechanical resonator with the frequency 
ujm- The parameter g describes the coupling strength between 
the cavity field and the mechanical resonator. The parame- 
ter n is the coupling strength between the cavity field and the 
external driving field with the frequency ujd and the phase (pd- 
If an unitary transform U = cx-p[ga^a{b'^ — b)/ujm] is ap- 
plied to Eq. (HJ, then the Hamiltonian in Eq. ^ with 17 = is 
diagonalized to hLUmb''b + hoja'^a — {hg'^a^'^a? jujm), where 
uj = ujc — g^ jijjra- That is, the mechanical resonator is 
decoupled from the cavity field when there is no driving 
field. Moreover, the nonlinear photon-photon interaction term 
hg^a^'^c? jujra induces the photon blockade |19] in the op- 
tomechanical systems with the strong coupling strength g and 
low dissipation of the cavity field. In this case, the driving 
field couples only two lowest energy levels |0) and |1) of the 
cavity field, and Eq. ([TJ with the transform \J can be reduced 
to the effective one 

iftw = h^o, + h^.^b'^b + h |n(t)(7+e''(^'-'') + h.c.} , (2) 

under the two-level approximation for the cavity field, where 
= fiexp \—i(ujdt + (/"d)]. Here, we redefine the photon 
operators via the Pauli operator cr^ = |1)(1| — |0)(0| and the 
ladder operator (T+ = |1)(0| ( cr_ = |0)(1|). Hereafter, we 
use |e) and \g) to denote the single-photon excited state |1) 
and the vacuum (ground) state |0) of the cavity field, respec- 
tively, i.e., |1) = |e) and |0) = \g). However, the states \k) 
with A; = 1, 2, • • • , denote the phonon number states of the 
mechanical resonator after the transform U . The parameter 
77 = g /uj,n is equivalent to the Lamb-Dicke parameter in the 
system of trapped cold ion. The third term in E g. (pl l, i.e., the 
driving-induced term, can be further written as 13911 

gt,h-^»(Oe"^^+I] .,„ +h.c., (3) 




FIG. 2: (Color online) Schematic diagram for the information leak- 
age to the third level due to the small arharmonicity. Here, two low- 
est horizontal lines in each column linked by the vertical black line 
with two arrows denote two-level approximation with the carrier pro- 
cess. The gray vertical line with the arrow pointed to the top line in 
each column simply denotes the information leakage in the carrier 
process. However, each red (dark red) slanted line with two arrows 
pointed to two black lines in different columns denotes the red side- 
band excitation (information leakage in the red sideband excitation). 



which corresponds to the carrier, red- and blue-sideband exci- 
tations in different conditions with the language of trap ions. 

As schematically shown in Fig.[Tl the carrier process links 
the transitions between |n, e) and \n,g) under the condi- 
tion ujd — oj, with an effective Rabi frequency i7„ — 
ncM-v'mj:;^,i-iyr^^^Ci/jl and Ci = nl/jiin - 
therefore, no phonon absorption or emission occurs, 
the external field only flips the photon states. The k- 
phonon red sideband excitation process links the transi- 
tions between \n,e) and \n + k,g) under the condition 
LUd ~ Ol! — kujm, with an effective Rabi frequency ftn.k ~ 
ncM-v'mvWi^ + k)\/n\Y';=^{-iyv^'CU{] + fc)!, 
i.e., the cavity field transits from the ground (excited) state to 
the excited (ground) state by absorbing (emitting) k phonons 
assisted by the external field. The fc-phonon blue sideband 
excitation links transition between |n, and \n + fc, e) un- 
der the condition a;^ = a; + fcwm, with an effective Rabi 
frequency Vln.k, i-c, the cavity field transits from the ground 
(excited) state to the excited (ground) state by emitting (ab- 
sorbing) k phonons with the help of the external field. Un- 
der the condition rf^n + 1 = g\/n + l/cj„i <C 1 with the 
average phonon number n of the mechanical vibration, we 
have cxp [77(6^ — b)] w 1 + 77(6^ — b). In this case, only a 
single-phonon transition occurs assisted by the driving field 
with fc = 1 in Fig.[T] The dynamical evolutions governed by 
the carrier, red and blue sideband precesses can be described 



by t/o=(t°), Ul and ^{t^), respectively ||39|], here the 



superscripts (e.g., r) denote different processes (the red side- 
band excitation process). 

Synthesizing phonon states. — We have reduced the Hamil- 
tonian in Eq. ^ of the driven optomechanical system to that 
in Eq. (|3]l, which is similar to one of trapped ions [8, SL\. 
Thus arbitrary superpositions — X]fc=o Ck\k) of different 
phonon states |fc) can be generated by using similar method as 
in trap ions [8,9]. As an example, below we show the detailed 
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procedure for generating the phonon state | k) . 



We assume that both the mechanical resonator and the cav- 
ity field are initially in the ground state |0) and the the vac- 
uum (ground) state \g), respectively. A single-phonon state 
|1) can be generated by two steps: (i) the cavity field is reso- 
nantly excited from the ground state \g) to the excited state 
|e) by the driving field with frequency matching condition 
= w and time interval — 7r/(2r2o,o)- That is, the 
optomechanical system is prepared in the state |e. 0) via the 
carrier process, (ii) The frequency of the driving field is set 
to = — ujm such that the red sideband excitation is im- 
posed by a time interval t\ = 7r/(2f2o,i), then the system is 
in the state [.g, 1), which means that the single-phonon state 
|1) is generated and the cavity field is in its ground state \g). 
Using the similar steps as for generating phonon state 1 1), any 
phonon state |fc) can be generated. However, we note that 
there are two ways to generate the state |fc) due to the ex- 
perimental constrain of the parameter 77. Under the condition 
ry <C 1, the generation of the state \k) requires /c-step carrier 
processes and fc-step red sideband excitation processes |8i], 
and the total time is X]i=o k/(2^^i,o) + 7r/(2rii,i)]. Out of the 
regime 7/^1, the state can be conveniently generated [9] 
by one step using the blue sideband excitation with the fre- 
quency matching condition oj^ = + kuj„i and the evolution 



time is = 7r/(257o.fc); or by two steps using one carrier 
(ujd — ^) and one red sideband excitation (w^ = lu ~ ktOm) 
processes for uj > kujm with total time t — + t'[, where the 
times for the carrier and red sideband excitation processes are 
<g = 7r/(2i7o.o) and = 7r/(2i7o,fc), respectively. It should 
be noted that the phonon state cannot be generated using 
one carrier and one red sideband excitation for u < kuj„i. 
More general discussions on synthesizing arbitrary superposi- 
tions of different phonon states |fc) can be found in sup- 
plementary materials 13911 . 



Leakage due to small anharmonicity of cavity field. — We 
know that the anharmonicity of the cavity field induced by 
the radiation pressure is not very large because the optome- 
chanical interaction is usually not strong. To study how the 
small anharmonicity affects the fidelity of generating nonclas- 
sical phonon states, we now study, as an example, the interac- 
tion between the phonon and the three-level photon system, 
formed by the three lowest energy levels |0) = \g), |1) = |e) 
and 1 2) = |e') of the cavity field. The transition frequency be- 
tween the states |e) and |e') is assumed as w + (5. The param- 
eter 6 characterizes the anharmoncity of the energy levels for 
the cavity field. The harmonic and the two-level model can be 
recovered when (5 = and 6 = 00, respectively. In optome- 
chanical systems, the anharmonicity is (5 = —2g'^/uj„i, which 
is a negative number, i.e., the transition frequency between the 
states |e) and |e') is smaller than that between the states \g) 
and |e). As schematically shown in Fig.|2l the Hamiltonian 
i/thi between the three-level photon system and the phonon 
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FIG. 3: (Color online) Fidelities for preparing states |2) in (a) and 
(|0) — |2})/\/2 in (b) using both numerical (black solid curve) and 
approximately analytical (red dash curve) approaches for rj = 0.1. 



mode can be written as 



iJthr = nwrab'^b + huj\e){e\ + h{2uj + 5)\e'){e'\ 



(4) 



[m [\e){g\^V2\e'){e 



H.c.j 



by projecting the cavity field operator in Eq. ([T]l with the uni- 
tary transform U to its three eigenstates jg), |e), |e'). Here, 
the parameter VL{t) in Eq. (HI is the same as that in Eq. (|2]). 

Let us first calculate the fidelity for preparing excite state 
|e) in the carrier process when the photon leakage to the third 
energy level is taken into account. In the carrier process for 
the time interval = 7r/(2ilo,o), the cavity field is prepared 
to its first excited state |e) from the ground state \g) under 
the two-level approximation in Eq. However, when the 
information leakage from the first to the second exited state 
is considered, the wavefunction of the cavity field should be 
written as |(^(t)) = Cg{t)\g) Ce{t)\e) + Ce' {t)\e') at the time 
t. Three coefficients Cg{t), Ce{t), and Ce'{t) can be obtained 
by solving the Schrodinger equation with given initial condi- 
tion (3^. Thus, under the condition fl <^ S, the fidelity of 
preparing the excited state |e) can be approximately given as 
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(5) 

If the ratio 77 is big enough such that the state \ k) can be gen- 
erated by a carrier process and a fcth red sideband excitation, 
then the phonon state \k) can be generated with the fidelity 
F as shown in Eq. (|5]l. However, if the ratio 77 is very small, 
then we need several steps to generate the state |fc). More- 
over, the preparation for arbitrary superpositions of different 
phonon states also needs several steps. Thus, the fidelity cal- 
culation, for the state preparation with several steps, becomes 
complicated when the information leakage is included. As 
examples, we study the leakage effect on generating the states 
|2) and (|0) — |2))/V2 by using the carrier and the single- 
phonon red sideband excitation processes. We have known 
that the generation of those states needs two carrier processes 
and two single-phonon red sideband excitations. The detailed 
calculations of the fidelity to achieve these two target states 
can be found in our supplementary materials [[3911 ■ Here we 
just give the approximated fidelities, Fi and F2 for preparing 
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TABLE I: Overview of the Lamb-Dicke parameters reached in microwave and optical cavity optomechanical experiments, where Uc and 7c 
are the frequency and the decay rate of the cavity, ujrn and 7™ are the frequency and the decay rate of the mechanical resonator, g is the 
optomechanical coupling constant, and r\ = gj i^m is the Lamb-Dicke parameter. 



System 


LJc/27r (Hz) 


7<:/27r (Hz) 


Qc = ^.h. 


uj„^/2-n (Hz) 


7„i/27r (Hz) 


Qm = cj™/7™ 


g/2TT (Hz) 


V = g/t^m 


Microwave cavity ^6] 


7.47 G 


170 K 


3.77 X lO" 


10.69 M 


30 


3.3 X 10^ 


226 


2.11 X 10"^ 


Membrane [32] 


282 T 


0.32 M 


8.81 X 10** 


134 K 


0.12 


1.12 X lO** 


2.68 


2 X 10"^ 


F-P cavity [3J1 


282 T 


215 K 


1.31 X 10^ 


947 K 


140 


6.76 X 10^ 


27 


2.85 X 10"" 


Toroidal microcavity [341 


385 T 


7.1 M 


5.42 X 10^ 


78 M 


10 K 


7.8 X 10^ 


3.4 K 


4.36 X 10"^ 


Optomeclianical crystals ^5] 


195 T 


500 M 


3.9 X 10^ 


3.68 G 


35 K 


1.05 X 10^ 


910K 


2.47 X lO""" 


Zipper cavity [36] 


194 T 


6G 


3.23 X 10* 


7.9 M 


98.75 K 


80 


599 K 


7.58 X 10"^ 


Double-wheel microcavity [37] 


190 T 


10 G 


1.9 X 10'' 


8.05 M 


2.01 M 


4 


732 K 


9.09 X 10"^ 


BEC [38J 


385 T 


1.3 M 


2.96 X 10** 


15.1 K 






0.39 M 


25.828 



States 1 2) and (|0) - |2))/V2, as 



F2 



frg if a 



f eg fee H~ frg {fee) 



(6) 
(7) 



under the condition ^ = 2g'^/ujm- Where /,.g = (1 — 



2ri^a)sm[Tr{4:-3rj'^a)/8], fee = (1 - 1.5a) sin 61, and /^g = 
(1 - 2a) sin 61 with a = n^/S^ and 61 = 7r(2 - a)/4. In 
Figs. [3a) and (b), Fi and F2 are numerically calculated and 
also compared with the approximated solutions in Eqs. ^ and 
(|7J. We can find that the fidelities tend to one when the rate 
is bigger than 20, moreover, the fidelities are bigger than 
0.9 v/henS/n > 10. 



Experimental requirements. — Let us now discuss the ex- 
perimental requirements for our proposal, (i) Similar to the 
ground-state cooling of the optomechanical system, the gen- 
eration of arbitrary superpositions of phonon states relies on 
the sideband excitations. This means that the frequency aj,„ of 
the mechanical resonator and the decay rate 7c of the cavity 
field have to satisfy the condition > 7c- (ii) Our proposal 
should work at the single-photon strong coupling regime as 
for the photon blockade [I9i, 2M in optomechanical systems, 
thus the coupling strength g and the frequency Wm of the me- 
chanical resonator should be larger than the decay rates jc 
and 7m of the mechanical resonator and the cavity field, i.e., 
'^mi 9 ^ Ic, 7m- Moreover, the nonlinear photon-photon 
interaction strength juj^ should be bigger than the decay 
rate 7c of the cavity field, i.e., (g^/cj„i) > 7c, such that the 
single-photon excitation or photon blockade can be guaran- 
teed and the two-level approximation can be applied, (iii) 
Negligible information leakage requires that the strength O 
of the classical driving field should be smaller than the an- 
harmoncity 2g'^ /tUm in the carrier process, (iv) As in trapped 
ion systems [gj, the big Lamb-Dicke parameter ry — gjuJm 
corresponds to the fast preparation of the multi-phonon states. 
Thus, (iii) and (iv) tell us that the big 77 is more desirable for 
our proposal. 

Based on the above requirements, we can estimate the ex- 
perimental parameters for realizing our proposal. For in- 
stance, if we take the moderate parameters 7c/27r = 0.1 MHz 
and 7,„/27r = 0.01 MHz, 5/27r = 10 MHz and w„i/27r = 20 
MHz, then phonon states |2) and (|0) — 1 2)) /a/2 can be per- 
fectly prepared when the strength of the driving field is taken 
as f7/27r = 0.5 MHz and 1 MHz, respectively. In Table U 



we have summarized experimental parameters for different 
setups of optomechanical systems. We find that the promis- 
ing candidates for realizing our proposal are the optomechan- 
ical crystal devices, superconducting circuits, and the ultrcold 
atoms in optical resonators with the current state of the art. 
Because these systems may have big ratio 77 and also satisfy 
all above requirements. Moreover, the condition 77 > 1 has 
been applied to generate superpositions of a mirror and pho- 
tons 121.1 in standard optomechanical systems. Thus our pro- 
posal should be feasible for experiments. 

Conclusions. — In summary, we have proposed to synthe- 
size arbitrary phonon states in optomechanical systems by us- 
ing the sideband excitations. Similar to the photon blockade 
1 19, [24], our proposal relies on the single-photon strong cou- 
pling condition such that the two-level approximation for the 
cavity field can be made in the optomechanical systems. We 
have shown that the parameters for our proposal are within the 
reach of the current technology. Our proposal opens up a way 
to engineer the arbitrary nonclassical phonon state on chip. 

YXL is supported by the National Natural Science Foun- 
dation of China under Nos. 10975080 and 61025022. JZ is 
supported by the NSFC under Nos. 61 174084, 61 134008. 
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